Manuscript for Orbital2001 



Mechanisms for Non- Trivial Magnetization Plateaux of an S =1 

Frustrated Spin Ladder 



O 

o 

m 



u 



X3 
O 

o 



> 
in 
cn 
o 
o\ 
o 



o 



X 
J3 



Kiyomi Okamoto, Nobuhisa Okazaki^ and Tom Sakai^ 

Department of Physics, Tokyo Institute of Technology, Oh-okayama, Meguro-ku, Tokyo 152-8551, Japan 
^ Kyoto Miyama High School, Murashita, Sasari, Miyama-cho, Kitakuwata-gun, Kyoto 601-0705, Japan 
^ Tokyo Metropolitan Institute of Technology, Asahigaoka, Hino, Tokyo 191-0065, Japan 

We investigate the non-trivial magnetization plateau at 1/4 of the saturation magnetization 
of 5 = 1 spin ladder, especially with reference to recent experimental results on a new or- 
ganic tetraradical 3,3',5,5'-tetrakis(A''-teri-butylaminoxyl)biphenyl, abbreviated as BIP-TENO. 
We propose three mechanisms for the formation of the plateau; the Neel mechanism, the dimer 
mechanism and the spin-Peierls mechanism. We also discuss the effect of four-spin exchange 
interactions. 
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§1. Introduction 

A new organic tetraradical, 3,3',5,5'-tetrakis(jV-tert- 
butylaminoxyl)biphenyl, abbreviated as BIP-TENO, has 
been synthesized very recently.B' This material can be 
regarded as an S = 1 antiferromagnetic two-leg spin lad- 
der. Goto et ala measured the magnetization curve of 
BIP-TENO at low temperatures in high magnetic fields 
up to about 50 T. They found M ~ up to about 
Hci — 10 T, suggesting the existence of the nonmag- 
netic ground state, which was consistent with the result 
of susceptibility measurement. Another remarkable fea- 
ture of their magnetization curve is the existence of the 
plateau at M = Ms/4 above Hc2 — 45 T, where Mg is 
the saturation magnetization. Unfortunately, the end of 
the Ms/4 plateau is unknown because of the limitation 
of the strength of their magnetic field. 
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Fig. 1. Simplified model of BIP-TENO. Open circles represent 
S = I spins. 

A simplified modefl) of BIP-TENO is shown in Fig. 1, 
although the real material is more complicated. Both 
coupling constants are thought to be antiferromagnetic. 
In experimental estimation, it seems Ji. 

Considering the necessary condition for the magne- 
tization plateau,EI' the spontaneous symmetry breaking 
along the leg-direction is necessary for the formation of 
the plateau. By use of the model of Fig. 1, we have al- 
ready proposed two mechanisms for thajplateau; the Neel 
mechanism and the dimer meclianism.Q|j^ 

In this paper, we briefly review the Neel and dimer 



mechanism, and also propose a new possible mechanism 
— the spin-Peierls mechanism. We also discuss the effect 
of four-spin exchange interactions. 

§2. Degenerate Perturbation Theory 

The Hamiltonian of Fig. 1 is expressed as 

L L 
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(2.1) 



where S denotes the spin-1 operator, j the rung number 
and 1 — 1,2 the leg number. The last term is the Zeeman 
energy in the magnetic field H. When J± = 0, our 
model is reduced to that of two independent conventional 
5=1 chains. In this case there appears no plateau in the 
magnetization curve except at M = 0, which corresponds 
to the Haldane gap. 

In the opposite limit J± ^ Ji we can use the degener- 
ate perturbation theory. When Ji = 0, all the rung spin 
pairs are mutually independent. Thus, at M = Ms/4, 
half of the rung spin pairs are in the state 



V'(o,o) = 
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and the remaining half pairs are in the state 



V'(l,l) = 



V2\ 



(2.2) 
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where tpi^tot, •S'tot) the wave function of with the quan- 



tum numbers S'tot and St, 



These two wave functions 



have the lowest energies in the subspace of S^^^ = and 
S'tot = respectively. The Ms/4 state is highly degen- 
erate as far as Ji = 0, because there is no restriction 
for the configurations of these two states. This degener- 
acy is lifted up by the introduction of Ji . To investigate 
the effect of Ji, we introduce the pseudospin T with 
T = 1/2. The I tr) and | i)-) states of the T spin corre- 
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spond to V'll, 1) and -0(0, 0), respectively. Neglecting the 
other seven states for the rung spin pairs, the effective 
Hamiltonian can be written as, through straightforward 
calculations 



Hcff — 



E 
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in the lowest order of Ji, where 



Jl = (8/3)Ji, 



J^^ - (l/2)Ji 



(2.4) 



(2.5) 



Thus, the Afs/4 plateau problem of the original model 
in a magnetic field is mapped onto the M = prob- 
lem of the T = 1/2 antiferromagnetic XX Z spin chain 
with the nearest-neighbor (NN) interaction in the ab- 
sence of magnetic field. As is well known, depending on 
whether Aoff = J^g/J^g < 1 or AoS > 1, the ground 



state of eq. (2.4) is either the spin-fluid state (gapless) 
or the Neel state (gapful), which correspond to the no- 
plateau state or the plateau state of the original model 
at Ms/4, respectively. Since Acg = 3/f6 in the present 



case, no Ms/4 plateau is expected as far as eq.(2.4) can 



be applied. As already stated, neither is there an Ms/ 4 
plateau when J± = 0. Thus, the Afs/4 plateau will not 
appear in the entire region of Ji / J± for the model of 
Fig. 1. We have checked this conclusion by numerical 
methods. 
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Fig. 2. (a) Physical picture of Neel mechanism plateau. Open and 
shaded ellipses represent ip{0, 0) and ip{l, 1) states of the rung 
pairs (I JJ.) and | in the T-picture), respectively, (b) Physical 
picture of dimer mechanism plateau. 

If we add frustrated spin interactions 
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to the Hamiltonian (2.1), the coupling constants in the 
effective Hamiltonian are modified to 



^eff =(8/3)(Jl- J2), 



J,^ff = (l/2)(Ji-|- J2), (2.7) 



and also there appear the next-nearest-neighbor interac- 
tions 
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Thus our effective Hamiltonian is of the generalized ver- 
sion of the T = 1/2 XX Z chain with the next-nearest- 
neighbor interactions. In this case there is a possibility 
of the Neel state and the dimer state of the T-system.t3' 
Although both states are gapped states (plateau states 
in the 5-picture) , the physical picture of the plateau for- 
mation is quite different with each other. 

We have already shownlf that the condition for the 
realization of the Neel plateau is J2/J1 > 13/19 when 
J3 = 0. Also the dimer plateau condition is J3/J1 > 0.31 
when J2 = 0. 

§3. Spin-Peierls Mechanism 

Another possible mechanism is the spin-Peierls mech- 
anism. As is well known, usual spin-Peierls transition 
is observed in 5* = 1/2 isotropic Heisenberg chain at 
M = 0. The uniform nearest-neighbor spin interac- 
tion JSj ■ Sj+i is modified into the alternating form 
J(l -I- Sy Sj ■ Sj+i due to the lattice distortion, where 6 
measures the magnitude of the bond-alternation. Since 
the dimer correlation behaves as in the long distance 
limit for the uniform case (quasi long-range order), the 
uniform chain is very sensitive to the bond alternation. 
In fact, the energy gain of the spin system due to the 
bond-alternation behaves asO S"' with a = 4/3, where 
a is the energy gain exponent.. Here we have omitted 
the logarithmic correction.!]^ Because the energy loss of 
the lattice system is proportional to 5^, there occurs the 
bond-alternation induced by spontaneous lattice distor- 
tion. Two spins on a stronger bond form a singlet pair 
effectively. Thus there is a finite energy gap in the ex- 
citation spectrum, which corresponds to the formation 
energy of an effective singlet pair. 
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Fig. 3. Physical picture of the spin-Peierls mechanism. The dis- 
tance between neighboring two spins on a leg changes alternat- 
ingly. This is the origin of the bond-alternation. 

Since the original Hamiltonian can be mapped onto the 
r = 1/2 XX Z chain, a similar spin-Peierls scenario is 
also possible for the present model. The physical picture 
of the spin-Peierls mechanism is shown in Fig. 3. The 
energy gap in the excitation spectrum of the T system 
corresponds to the magnetization plateau at Ms/4 of the 
original system. Of course, the value of the energy gain 
exponent a should be modified from 4/3, because our 
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effective Hamiltonian is no longer isotropic (now Acff = 
3/16 7^ 1). The Acff dependence of a isffl) 

4 2 

V = T-T^TT-fz-nr- (3.1) 



4-r]' 
which results in 



1 + (2/7r)cos-i Aeff 



1.80 < 2 



(3.2) 



for our case. Thus the spin-Peierls mechanism is possible 
from the consideration of the energy. Our preliminary 
numerical calculation based on the 5 = 1 model shows 



that the a is very close to 1.8 as conjectured in eq. (3.2) 



§4. Effect of Four-Spin Cyclic Exchange Inter- 
actions 
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Fig. 4. Four-spin cyclic exchange interactions. 

Let US consider the effect of four-spin cyclic exchange 
interactions shown in Fig. 4. In the 5* = 1/2 two- leg 
ladder, this interactiorL-can brings about the non-trivial 
magnetization plateau.EiP The four-spin cyclic exchange 
interaction works as 
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where a, (3, 7, S are the related four spins. If we add this 
interaction with magnitude J4 to the Hamiltonian (2.1), 
the effective coupling constants in eq. (2.5) are modified 
to 



„y 4J4 3Jl + 2J4 

'^off = Z ' "'off 



3 ^ 6 - ^'-'^ 

Here we do not take the J2 and J3 interactions into 
consideration for simplicity. It is easy to see that the 



Neel state in the T-system is not realized for any choice 
of J4. Thus the four-spin cyclic exchange cannot bring 
about the plateau in our problem. This situation ia_nuite 
different from that in the 5 = 1/2 two-leg ladder .I12P 

§5. Summary 

We have proposed three mechanisms for the plateau 
of 5 = 1 ladder, BIP-TENO. These mechanisms are di- 
rectly distinguished by experimental methods in prin- 
ciple, although they might be difficult to realize. In 
the Neel plateau state, the expectation values of 5^ of 
spins along a leg are 1/2, 0, 1/2, 0, 1/2, • • •, although this 
"complete Neel configuration" is somewhat smeared by 
the quantum fluctuations. On the other hand, in the 
dimer plateau state, they are 1/4, 1/4, 1/4, 1/4, 1/4, • • •. 
In the spin-Peierls plateau, the spontaneous lattice dis- 
tortion can be observed, although the spin configurations 
is the same as those in the dimer plateau state. 

Because such experiments seem to be difficult in high 
magnetic fields, we have to surmise the nature of the 
plateau state from the behaviors of the magnetization 
curve and the magnetic susceptibility, by comparing 
them with the analytical and numerical results. The in- 
vestigation in this line is now in progress. 



[1] K. Katoh, Y. Hosokoshi, K. Inoue and T. Goto: J. Phys Soc. 

Jpn 69 (2000) 1008. 
[2] T. Goto, M. I. Bartashevich, Y. Hosokoshi, K. Kato and K. 

Inoue: Physica B 294 & 295 (2001) 43. 
[3] M. Oshikawa, M. Yamanaka and I. Affleck: Phys. Rev. Lett. 

78 (1997) 1984. 

[4] K. Okamoto, N. Okazaki and T. Sakai: J. Phys. Soc. Jpn. 70 
(2001) 636D 



[5] N. Okazaki, K. Okamoto and T. Sakai: cond-mat/0108527 



[6] M. C. Cross and D. S. Fisher: Phys. Rev. B 19 (1979) 402. 
[7] J. L. Black and V. J. Emery: Phys. Rev. B 23 (1981) 429. 
[8] T. Nakano and H. Fukuyama: J. Phys. Soc. Jpn. 50 (1981) 
2489. 

[9] H. Nishimori and K. Okamoto 
1132. 

[10] K. Nomura and K. Okamoto: 
(1994) 5773. 

[11] A. Nakasu, K. Totsuka, Y. Hasegawa, K. Okamoto and T, 
Sakai: J. Phys.: Condens. Matter 13 (2001) 7421. 



J. Phys. Soc. Jpn. 56 (1987) 
J. Phys. A: Math. Gen. 27 



